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Abstract 

We consider quantum systems with a chaotic classical limit that are perturbed by a point-like 
scatterer. The spectral form factor K(t) for these systems is evaluated semiclassically in terms 
of periodic and diffractive orbits. It is shown for order r 2 and r 3 that off-diagonal contributions 
to the form factor which involve diffractive orbits cancel exactly the diagonal contributions from 
>► . diffractive orbits, implying that the perturbation by the scatterer does not change the spectral 

•i-H . 

^ ■ statistic. We further show that parametric spectral statistics for these systems are universal for 
c$ • small changes of the strength of the scatterer. 
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1 Introduction 



Semiclassical theories for spectral statistics have been developed 0, §, || to find an explanation for 
the observed universality in energy spectra of quantum systems with a chaotic classical limit, the 
agreement of correlations in energy spectra with those between eigenvalues of random matrices 
[f|]. They are based on semiclassical trace formulas that approximate the density of states in terms 
of classical trajectories ||. It has been shown by these theories that in the asymptotic limit of 
long-range correlations two-point correlation functions do coincide with those of random matrix 
theory 0, U . These results are based on mean properties of periodic orbits |J . To go beyond the 
leading asymptotic term requires information about correlations between periodic orbits which 
are presently not available |J. 

One of the expectations, on basis of the random matrix hypothesis |7j , is that a perturbation of 
a chaotic system should not change the statistical distribution of the energy levels of the quantum 
system, if it does not change the chaotic nature of the classical motion. In the present article we 
investigate, on the level of the semiclassical approximation, whether the perturbation by a point- 
like scatterer has this property. One argument in favour of this invariance is that the semiclassical 
approximation for the density of states is not changed in leading order of h for this perturbation. 
The influence of the scatterer is described semiclassically by a certain class of trajectories, so-called 
diffr active orbits that start from the scatterer and return to it. They contribute to the density of 
states in higher order of h than the leading order contribution from periodic orbits. 

The present article is motivated by the observation in || that a scatterer could neverthe- 
less have an influence on spectral statistics. When spectral correlation functions are calculated 
by using mean properties of diffractive orbits, the so-called diagonal approximation, they show 
modifications which, in general, do not vanish in the semiclassical limit {K — > 0). In order that 
this does not lead to deviations from random matrix statistics, these terms have to be cancelled 
by off-diagonal terms which contain information about correlations between different trajectories. 
As remarked above, the calculation of correlations between trajectories is an unsolved problem 
in general systems. For the diffractive orbits that describe the influence of a scatterer, however, 
off-diagonal terms can be calculated explicitly. This is done in the following sections. The results 
show that diagonal and off-diagonal terms indeed cancel each other. Furthermore, the results can 
be used to investigate parametric spectral correlations, i. e. correlations between spectra of the 
system for different parameter values, where the parameter is the strength of the scatterer. It is 
shown that the parametric spectral correlations are universal for small changes of the parameter. 

2 The spectral form factor 

The perturbation by a point-like scatterer is represented, formally, by a delta-potential H = 
H + X5(r — r ), where A and r are the strength and position of the scatterer, respectively. In 
more than one dimension such a delta-potential is, however, not well defined. For example, it 
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Figure 1: Example of a double-diffractive orbit. 

leads to a divergent expression for the Green function. The problem can be regularised by the 
method of self-adjoint extensions, leading to a one-parameter family of Hamiltonians. A detailed 
monograph with references on the history and on applications of delta-like potentials is |J. We 
use in the following the property that the semiclassical approximation for the density of states 
has the same form as in the geometrical theory of diffraction 11], [12[ (see also |14|, [15| for 
applications on spectral statistics). 

We consider chaotic systems whose Hamiltonian is given in terms of a scalar and a vector 
potential. Billiard systems can be included in this description by letting the scalar potential be 
infinite outside the billiard region. The statistical distribution of the energy levels is investigated 
by semiclassically approximating the spectral form factor. We restrict to two-dimensional systems 
in order to keep the notation simple, but analogous calculations can be performed in higher 
dimensions. 

The spectral form factor is defined as Fourier transform of the spectral two-point correlation 
function 

K(r) = j°° Jj?L ^ osc (e + |) d osc (E - |) ) e exp (2<Ki W d{E)) . (1) 

The function d osc (E) = d(E) — d(E) is the oscillatory part of the density of states, and d(E) is 
the smooth part which is given in two dimensions by d(E) ~ T l (E)( y 2iih)~ 2 , E — > oo, where £(-£7) 
is the volume of the surface of constant energy in phase space. The statistics is evaluated by 
averaging over an energy interval that is small in comparison to E but contains a large number 
of energy levels. 

The semiclassical approximation for K (r) is obtained by inserting into (|1|) the approximation 
for the oscillatory part of the density of states 



dosc(E) w ^: Re ^ A, exp (^(E)j . 



(2) 
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In systems with a delta-like potential the sum in @ runs over all periodic orbits 0, and further 
over all diffractive orbits that start from the scatterer and return to it an arbitrary number of 



times n [Hi O, O, [HJ. An example for a double-diffractive orbit (n = 2) in a billiard system is 



shown in figure [l]. For n-fold diffractive orbits the amplitude has an ^-dependence of h n ^ 2 , and 
Sj denotes the action of an orbit. 

With (|2]) one obtains the following approximation for the spectral form factor 

where T = 2Trhd(E)r, and T 7 is the period of an orbit. For small values of r one can evaluate the 
double sum in (^) in the diagonal approximation ||. One obtains in this way from the periodic 
orbits the correct random matrix result K{r) ~ |r, r — > 0, where (3 = 1 or 2 for systems with or 
without time-reversal symmetry, respectively. 

The diagonal contributions from diffractive orbits to the form factor have been calculated in 
H . The result for n-fold diffractive orbits is 

( \ \V\ 2n T n+l 



where T> is the diffraction coefficient for the diffraction on the singularity of the potential [§, |l| 
It can be parameterised in the following form 

V = — T~r~\ ■ (5) 

if -7 -log (&) 

Here k = ^/2m(E — V(ro))/h, Vq is the position of the scatterer, a is a parameter describing the 
strength of the potential, and 7 is Euler's constant. In order that the terms (f|) do not lead to a 
deviation from random matrix statistics they have to be cancelled by off-diagonal terms involving 
diffractive orbits. By calculating off-diagonal terms for order r 2 and r 3 we show in the following 
that such a cancellation does indeed occur. 

We note that the diffraction coefficient satisfies the identities 

|£f = -4Im£> , |£>| 4 = 8(|£f -ReV 2 ) , (6) 

that will be used in the following. The first of these relations expresses the conservation of 
probability, and the second is a consequence of the first one. 



3 First-order correction 

The first-order correction to the diagonal approximation for the form factor arises from off-diagonal 
terms in @ between periodic orbits and single-diffractive orbits. In leading order, these orbits 



4 



Figure 2: A diffractive orbit (full) that is almost periodic and a nearby periodic orbit (dashed). 
The local coordinate system is oriented along the final direction of the diffractive orbit. 





are only correlated if the diffractive orbit follows the periodic orbit very closely. This happens, if 
the diffractive orbit is almost periodic, i. e. if the final momentum is almost identical to the initial 
momentum. An example is shown in figure |2|. 

The periodic orbit can be described by linearising the motion around the diffractive orbit. The 
condition that a trajectory in the vicinity of the diffractive orbit is periodic leads to the following 
equation 

/ X \ I X \ 

(7) 

Here M is the stability matrix of the diffractive orbit, e is the angle between the initial and final 
direction of the diffractive orbit, 7 is the angle between the direction along the periodic orbit 
and the final direction of the diffractive orbit, and 5 is the spacial distance between periodic and 
diffractive orbit (see figure H). The quantity p v is defined as mass times velocity, p v = mv. The 
index v is used in order to distinguish it from the canonical momentum p in systems with magnetic 
field. The stability matrix for the motion in a magnetic field is discussed in the appendix. 

In the linear approximation the difference in actions is obtained by expanding the action up 
to second order 

AS(E) = S po (E) - S ldo (E) « -pjs + ±(Ap f - A Pl )5 = -±5ep v , (8) 

where Apf and Api are the differences between the initial and final momenta of the periodic orbit 
and the diffractive orbit, respectively. The solution of the linear equation (|7]) yields the following 
relation between 5 and s 

M 12 

= — — — ep v , (9) 

Tr M — 2 w 

so that AS(E) depends quadratically on e. 
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With this approximation the off-diagonal terms are calculated. The amplitude of the diffractive 
orbit is given by 



T^V I 2lth \ ,7T .37T 



A ™=^imM^v-^-^)' (10) 

where T 7 is the time along the orbit, M 7 is its stability matrix, and v 1 is the number of conjugate 
points along the orbit. For the periodic orbit the corresponding amplitude is 

NL° = . T \ exp{-4// 7 | . (11) 
7 ^|TrM 7 -2| I 2^ 7 J 

The stability matrix is the same for both orbits in leading order, but the Maslov index of the 
periodic orbit can differ from the number of conjugate points z/ 7 by 1 || 



/i 7 = ^ 7 + -(l-K 7 ) , K 7 = sign^ TrM ^_ 2 j . (12) 

In the following we sum over all diffractive orbits that are almost periodic. This is done by 
applying first the sum rule for diffractive orbits for which the angle difference between initial and 
final direction has a fixed value s, and then integrating over the angle e. The sum rule is given by 

i 

V (£) - \ ,^(T-T 7 )^^4, (13) 
^ |(M 7 ) 12 | 1 ,} E(E) 
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where S(-E) is the volume of the energy shell. (It is implied in (|H|) that the left-hand side is 
smoothed over small intervals of T and e in order to obtain a non-singular expression.) 

Finally, one has to determine the multiplicity factor of the contribution. First, each off-diagonal 
term in (^) has a corresponding complex conjugate partner. If the summation is carried out over 
only one of these terms one has to take twice the real part of the sum. Furthermore, the periodic 
orbit and the diffractive orbit both have multiplicities 2/3 -1 , but a particular constellation occurs 
2/3~ 1 times in the sum over e (for systems with time-reversal symmetry for e and — e), so the total 
multiplicity is g = 2/3~ 1 . 

Inserting the amplitudes (0) and (|TT|) , and the action difference (§) with (§) into (|3|) we obtain 

= di(I> jySf^Wr exp (-iAS,(B)) S(T - T 7 ) 

4 r°° (e) V \^2nhS(T — T~) ^Bk_,£ (2+): ) 

— RC / dc ;> , TrM 7 -2 1 4 { - Z+K 1> 

2nM(E)f3 ./.oo 4^ 4vr^ v /|(M 7 ) 12 (TrM 7 -2)| 



1 ^r^^ VymKT ~ Ti 



2irhd(E)(3 J.^ ^ 47rp w |(M 7 )i2 



a 2 

2h 



2T 2 /-oo VpvV ^h .'V 

lm / de 



2M{E)P 7_oo 
2r 2 

— ImD. (14) 
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The integration over e can be carried out from minus to plus infinity since the main contribution 
comes from the vicinity of e = 0. Furthermore, the following steps have been carried out. First the 
integration variable has been changed to make the exponent independent of the stability matrix. 



Then the sum rule (13 ) has been applied, assuming that the distribution of angles between initial 
and final momenta of a diffractive orbit is independent of the distribution of the elements of the 
stability matrix 



£ KM^^-^-E W^i S{T - T -' h £=£ VtVM^2' (15 » 

and finally the integration has been carried out. 

K^(t) is the leading order correction to the diagonal approximation for the form factor and 
it cancels exactly the diagonal contribution from single diffractive orbits ((f|) with n = 1). This 
can be seen by using (^) 

K?\r) + k2>{t) = |^r 2 + |r 2 ImV = . (16) 

It shows that the presence of a point-like scatterer does not modify the spectral form factor up to 
order r 2 in systems with a chaotic classical limit. 

In order to find the geometries of orbits which contribute to a given order in r it is helpful to 
count the orders of h. The m-th order off-diagonal correction to the form factor is a r m+1 -term 
with a coefficient that has to be h- independent. The prefactor of the double sum over orbits 
in (HD is of order frr 1 and the product of the amplitudes of a rii-fold and a 77-2-fold diffractive 
orbit is of order fj-( n i+ n 2)/ 2 5 w h ere periodic orbits are denoted here as 0-fold diffractive orbits. 
The conversion of time T m+1 into r m+1 gives an order h m+1 , which yields altogether an order of 
j^(2m-n 1 -n 2 )/2 _ Furthermore, every integration over a small parameter e gives an additional order 
fir 1 / 2 , if the action difference is quadratic in this parameter. As a consequence, 2m — ni — n 2 small 
parameters are necessary in order that the prefactor of 

r (m+i) is 

^-independent. For the first order 
correction in this section {m = ni = 1 and n 2 = 0) this estimate gives one small parameter e. 

4 Second-order corrections 

For the second-order corrections we consider orbits that return twice to the region in coordi- 
nate space from which they started. These orbits are close to double- diffractive orbits. Double- 
diffractive orbits have the semiclassical amplitude 

>2 do ^V 2 1 f .7T, . .37T 

■ml l\l D \ Tl 

^ ) 12 



where T 7 is the total time along the trajectory, L 7 and _R 7 are the stability matrices for the two 
loops ('left' and 'right'), and v 1j l and z^ 7 ,r are the number of conjugate points along the loops. 
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fllTf) is the amplitude for one particular sequence in which the loops are traversed. In systems 
without time-reversal symmetry the degeneracy of the trajectory is thus two, meaning that there 
is another trajectory with exactly the same semiclassical amplitude and action. This trajectory 
traverses first the second loop of 7, and then the first loop of 7. In systems with time-reversal 
symmetry the degeneracy is eight. 

The sum rule for double- diffract ive orbits is given by || 

if there are n restrictions to the four directions of the velocities at the point from which the 
trajectories start and to which they return. As will be seen in the following, it follows from 
this sum rule that the contributions are of order r 3 (there is a factor T from every semiclassical 
amplitude, and a factor T from the sum rule). 

There are several possibilities in which a double-diffractive orbit can have an action which is 
almost identical to the action of a single-diffractive or a periodic orbit. A necessary condition is 
that there is always at least one small relative angle between the different initial and final directions 
of the orbit at the scattering point. In order to find the relevant cases one has to consider all 
possibilities and take into account the ^-argument that was given at the end of section |3|. The 
result is that there are three relevant configurations for systems without time-reversal symmetry 
and five configurations for systems with time-reversal symmetry. They are discussed in the next 
sections. 



4.1 Correlations between double-diffractive and single-diffractive or- 
bits 

Correlations between double-diffractive and single-diffractive orbits exist if the double-diffractive 
orbit is almost single-diffractive. This occurs if the final velocity of one loop deviates by a small 
angle e from the initial velocity of the other loop. An example is shown in figure |^. There is only 
one small parameter here which agrees with the estimate 2m — n\ — n 2 for m = 2, m = 2 and 
n 2 = 1. 

The further calculations are done analogously to the last section. The motion in the vicinity 
of the double-diffractive orbit is linearised and one obtains in this approximation a condition for 
the neighbouring single-diffractive orbit 

' (")=*( / <«» 

yPvl2j \PvllJ \Pvl3j \Vv{l2-£)j 

The angles 71, 72, 73 and the distance S are shown in figure |3[ The difference in action is obtained 
by expanding the action up to second order 

AS(E) = S ldo (E) - S 2do (E) = ~6e Pv = %^ , (20) 



2 ru 2 M- 
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Figure 3: A double-diffractive orbit (full) for which one initial direction deviates by a small angle 
e from one final direction, and a nearby single diffractive orbit (dashed). 



where M = RL is the stability matrix of the single-diffractive orbit. The last step in (|20D follows 
from the solution of ( |19|) for 5. 

The number of conjugate points v 1 along the single-diffractive orbit can differ from the sum 
of the number of conjugate points along the two loops, v 7j l and f 7 ,R, by 1. The general condition 
for this is 

z/ 7 = Vlh + VlR + -(1 - (j 7 ) , a 7 = sign ^ — — j . (21) 

The single- and double-diffractive orbits have multiplicity 2/3 -1 and 8/?~ 2 , respectively, but each 
configuration occurs for 2/3~ 1 different values of e. Therefore the total multiplicity is g = 8f3~ 2 . 
The contribution to the form factor from all pairs of orbits is thus given by 

- dfe 2Re />£ w W°>* «p (-s AS ^ £ ») a ' T - T ^ 

16 r°° ( e) hT 2 V 2 V* V2nH5(T-T~) (l i)i2{ r i)i2 . 

= i^— -Re / deV ' e~ (^h 2 ^ {2+a -> 

2irhd(E)(3 2 J.^ ^ 647r 2 pVI( L 7)i2 TfQh* WM 

1 „ ^(s')h 2 T 2 V 2 V*5(T -T 7 ) 

Re^ — — ; ,„ N — ^-e 2 



2nhd(E)[3 2 ^ 2tt^|(L 7 ) 12 ( j R 7 ) 12 | 



r 



— lm(V 2 V*). (22) 
Here we slightly abbreviated the procedure of ([14]) and performed the integration directly. 

4.2 Correlations between double-diffractive orbits and periodic orbits 

In order that a double-diffractive orbit is close to a periodic orbit it has to be almost periodic. 
This means that the final direction of each loop has to be almost identical to the initial direction 
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Figure 4: A double- diffract ive orbit (full) for which the two initial directions deviate by small 
angles E\ and e 2 from the two final directions, and a nearby periodic orbit (dashed). 



of the other loop, so there are two small relative angles as is shown in figure f|. 

The linearisation of the motion in the vicinity of the double-diffract ive orbit leads to the 
following condition for the periodic orbit 

\ / x \ / x_ \ 

(23) 

The angles 7^ and distances <5j are defined analogously as before in terms of the local coordinate 
systems that are oriented along the two final directions of the diffractive orbit. 
The difference in actions is given by 







AS(E) = S po (E) - S 2do (E) 
= -y(5i£! + 8 2 e 2 ) 

p 2 v {RL) 12 e\ + {LR) 12 e\ + 2(L 12 + R 12 ) e x e 2 



2 Tr M - 2 



(24) 



where M = RL is the stability matrix for the periodic orbit. Equation (|24"D can be written in the 
terms of a symmetric matrix A such that AS(E) = —\v 2 vl^ij-^-ij e i e r The number of negative 
eigenvalues of A is given by the number of sign changes in the sequence of sub-determinants 
(1, An, det A). With A u = M 12 /(TrM - 2) and det A = L 12 R 12 /{TtM - 2) one finds that the 
two signs of the eigenvalues are given by 

. ( M 12 \ . ( L 12 R l2 \ 

k = sign — - , a = sign — — . (25) 

6 V TrM - 2 / \ M 12 J K J 

The Maslov index of the periodic orbit can now differ by zero, one or two from the sum of 
conjugate points along the left and right loop. The criterion for this is the combination of (O) 
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and (pi ) and has the form 



li = v L + v R + -{2-K-<r), (26) 

where k and a are given in (|25|). 

The multiplicities of double-diffractive and periodic orbits are 8/5~ 2 and 2/3 -1 , respectively, but 
a particular configuration of them occurs 4/? _1 times in the integral over the angles (for example, 
for systems without time reversal symmetry the two angles can be interchanged), so the total 
multiplicity is g = Af3~ 2 . 

After inserting the amplitudes (O) and flTTp and the action difference ( |2~3| ) into ([5|) we obtain 

K^\r) = ^py2Re J^de^ 1 '^ A^{A^)* exp (-^AS y (E^j *(T-T 7 ) 

^e 2 )hT 2 V 2 5{T - T 7 ) e# S A Wj -;f (4+^+^) 



/oo 
dei d£ 2 y~] 



2tt^(E)/? 2 ^ ^ 16^^/1(^7)12 W12 (TrM 7 -2) 



1 ^(e' 1 ,e' 2 )h 2 T 2 V 2 5(T -T, 



2M(E)^ ^ p4|(L 7 ) 12 (i? 7 ) 



12 



r 



3 



-Re(P 2 ). (27) 

4.3 Correlations between pairs of single-diffractive orbits 

For exactly the same kind of double-diffractive orbits as in the last section, there is one further 
type of correlation that has to be considered. It occurs because there are two possible ways 
in which the double-diffractive orbit can be deformed into a nearby single-diffractive orbit, and 
consequently, there are correlations between these single-diffractive orbits. 

The action difference between the two orbits can be obtained from the action difference between 
each of these orbits and the double-diffractive orbit (|20|) 

AS(E) = S^\E) - S^ 2 (E) = -U^ + U 2 e 2 p v = (^e? - , (28) 

where M = RL and iV = LR are the stability matrices of the two single-diffractive orbits. The 
number of conjugate points along the orbits are given by 

vi = v L + v R + -(1 - a x ) , V2 = u L + v R + -(1 - a 2 ) , 

ai = sign {-n^~ ) ' a2= Slgn vn^r ) ■ (29) 

The multiplicities of the two orbits is both 2f3~ l and, as before, the configuration occurs 4/3 -1 
times in the double integral over the angles, so the total multiplicity is g = (3~ l . 
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Figure 5: A double-diffractive orbit (dotted) for which the two initial directions deviate by small 
angles E\ and e 2 from the two final directions, and two nearby single-diffractive orbits (full and 
dashed) . 



Inserting the amplitude (p~Q|) with stability matrix M and N, respectively, and the action 
difference ( p8|) into (|3]) one obtains 

2 r°° (e 1 ,e 2 )2nhT2\V\ 2 5{T -T^e-^'W-^' 1 -^ 



/oo 
dei de 2 y^ j 
"°° 7 



27rW(J5))9 i-00 V 167r 2 pVI( M 7)i2 



1 „ ^i^)fr 2 T 2 |X?| 2 5(T-T 7 



2ttM(E)(3 ^ 2pi\(L 1 ) 12 (R r ) 12 \ 

V\ 2 . (30) 



r 3 

' '-— 12 



2/5 

The contributions ^ o) (r) ^g 6) (r) and K^ c) (r) are the only second-order off-diagonal cor- 
rections in systems without time-reversal symmetry. As will be shown in the following, these 
contributions cancel exactly the diagonal term K^\r) for (3 — 2. For systems with time-reversal 
symmetry there are further contributions. They arise from the possibility that one trajectory can 
follow one loop of another trajectory in the same direction, but the other loop in the time-reversed 
direction. We assume in the following that the relevant contributions come from orbits which are 
close to double diffractive orbits in coordinate space and we evaluate their contributions in the 
next two subsections. 
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Figure 6: A double-diffractive orbit (dotted) with one loop for which the initial direction and 
the reversed final direction deviate by small angles E\ and £2 from the final direction of the other 
loop. Furthermore, two nearby single-diffractive orbits (full and dashed), one of which traverses 
the second loop in the opposite direction (full). 



4.4 Correlations between pairs of single-diffractive orbits involving 
time-reversed loops 

The first possibility involves two single-diffractive orbits. These orbits follow closely a double- 
diffractive orbit and traverse one loop in the same direction and the other loop in the opposite 
direction. In order that this can occur there must be one loop which has a very small opening 
angle, and it must be almost aligned to the final direction of the other loop as shown in figure 

The action difference and the indices are given by the equations (pS|) and ( p9| ) but now with 
M = RL and N = R l L where R % is the stability matrix for the time-reversed second loop. In 
term of the elements of R the matrix R l is given by 

* = { f- R A ■ (3D 

\iV2l till J 

The multiplicity of each orbit is two, and the configuration occurs two times in the integral 
over the angles, so the total multiplicity is g = 2. Inserting the amplitude flIU|) with stability 
matrix M and N, respectively, and the action difference ( 28) into (|3|) results in 



^ )(T) = 2^b) /^^E^'V ' 2 ^ ' 1 )* exp (4AS 7 (£)) S(T-T 7 ) 
4 r°° ^(e 1 ,e 2 )27ThT^\V\ 2 5{T -T 1 )e-i As ^ E) ' i ^ (a ^ 1 ~ (7 ^ ) 



1 v^K,4) h 2 T 2 \V\ 2 S(T - T 7 



2M{E) Y ^1(^7)12(^7)12! 
r 3 \V\ 2 . (32) 
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Figure 7: A double-diffractive orbit (dotted) for which all initial and final directions lie almost 
in one line. Furthermore, a nearby periodic orbit (dashed), and a nearby single-diffractive orbit 
(full) which traverses the second loop in the opposite direction. 



4.5 Correlations between single-diffractive orbits and periodic orbits 
involving time-reversed loops 

The last relevant configuration occurs if all initial and final velocities of the double-diffractive orbit 
lie almost in one line as in figure [7]. Then there exist neighbouring single-diffractive and periodic 
orbits which follow one loop in the same direction and the other loop in the opposite direction. 

The action difference between the periodic orbit and the diffractive orbit is obtained from (p0|) 
and (H) 

AS(E) = S po (E) - S ldo (E) 

= - + 5 2 e 2 ) + -£S 3 e 3 

_P 2 V {RL) l2 e\ + {LR) 12 e\ + 2(L 12 + R 12 ) e x e 2 1 2 2 L 12 R 12 
"2 TrM-2 + 2 £3Pv N 12 ' {66) 

where M = RL and N = R l L are the stability matrices of the periodic and diffractive orbit, 
respectively, and the indices v and fi of the diffractive and the periodic orbit are 

v = v L + v R + ^(1 - a 2 ) , /i = v L + v R + ^(2 - k - a i) , 

. (L 12 R l2 \ . f M 12 \ . (L l2 R l2 \ 

o 2 = sign , k = sign , <Ti = sign . d4 

V N 12 J \TtM-2J ' & \ M 12 J K 1 

For each double-diffractive orbit like the one in figure |^ there are two periodic orbits of multi- 
plicity two and four single-diffractive orbits of multiplicity two in the vicinity, but only half of the 
possible pairs involve exactly one time reversed loop which makes a total of 16. Furthermore, the 
double-diffractive orbit occurs 8 times in the integral over the angles and the total multiplicity is 
thus g = 2. 
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Inserting the amplitudes fllOD and (|TT|) and the action difference (^) into ([|) yields 

Kfi\r) = ^py2Re J^de^^^ exp (-^AS 7 (E^j 8{T — T 7 ) 

4 r°o ^( ei , e2 , e3 )r2pv / 2^5(T-T 7 )e-i A ^^-^( 2 - CT ^ + ^ + ^ 1 ) 

= =- — -Re / deid£ 2 d£:3> — = 

2M{E) ^ 4vr^ v /|(M 7 ) 12 (TrM 7 -2)| 

1 ^(4,4,4) 47rn 2 T 2 Z> 5{T - T 7 ) _. f 

= 2M{E) G ^ piKL^R^l 6 

= 2r 3 Im£>. (35) 
The sum of all contributions can be written in the form 

K^\t) = Kf\r) + K%\t) + K%\t) + K%\t) + K™(r) + K%\t) 

T — (-|P| 4 + \V\ 2 ImD - Re£> 2 + -|D| 2 + (2 - (3)\V\ 2 + 2(2-/3) ImP 



/3 2 V8 2 
r 3 /l 

= — f -|P| 4 + \V\ 2 ImV - ReV 2 + \V 
= , (36) 

which can be seen by using (H). This shows that off-diagonal terms cancel the diagonal term 
also in this order. It implies that the form factor is determined by periodic orbits alone, because 
the different terms which involve diffractive orbits cancel each other. This might be true also for 
other point-like sources of diffraction like e.g. Aharonov-Bohm flux lines (see flT]]), although a 
quantitative analysis would require here the use of uniform approximations. 



5 Universality in parametric correlations 

The cancellation of off-diagonal and diagonal terms is conform with the expected universality of 
spectral statistics in chaotic systems. Universality is, however, not only expected in the properties 
of single systems, but also in the way in which system properties vary when a parameter of 
the system is changed. For example, random matrix theory makes predictions about correlations 
between densities of states for different parameter values. The semiclassical calculation of diagonal 
and off-diagonal terms allows to test this prediction for systems with a point-like scatterer where 
the parameter is the strength of the scatterer. 

In analogy to the spectral form factor, a parametric form factor can be defined as Fourier 
transform of the parametric two-point correlation function 

K(t, x) = J|y (ciosc (E + |, X + X - ) d osc (E - |, X - |) ) e exp (2m V rd(E)) . (37) 

Here x is the parameter difference between two systems. In order for this statistics to be universal 
the parameter has to be chosen in a particular way. The requirement is that the variance of the 
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velocities, the derivatives of the unfolded energies with respect to the parameter (de n /dx), is equal 
to unity. The unfolded energies are obtained from the quantum energies of the system by a scaling 
that leads to mean level distance of one. 

For random matrix ensembles, the parametric two-point correlation function was derived in 



18 1 in the context of disordered metallic systems. For the GUE-result, which we discuss first, the 



Fourier transform in (37) can be evaluated in a closed form. It results in 



{ sinh ( 27r x T ) exr,(-2ir 2 r 2 T) if r < 1 

27T 2 X 2 T exp^ZTT x T) it r < 1, 

sinh(27r 2 iE 2 T) / o 2 2 2\ -r ^ i 

27r 2 x 2 T ' exp(-2ir z x z T z ) if r > 1, 

and has for small values of r the expansion 

iT GUE (r, x)=r- 2ttVt 2 + 2vrW + . . . . (39) 

We examine in the following whether the perturbation by a point-like scatterer leads to univer- 
sal correlations. For large parameter differences x the parametric correlations for these systems 
cannot be expected to be universal. The treatment of a delta-scatterer by the method of self- 
adjoint extensions leads to a quantisation condition with the property that there is exactly one 
eigenvalue of the perturbed system within each pair of neighbouring eigenvalues of the unperturbed 
system. This puts a restriction to the movement of eigenvalues when the parameter is changed. 
For this reason, one can expect universal properties only for small parameter differences. 

We choose first a in ([5]) as parameter of the system. The two densities in (|37j) differ then 
only in the diffraction coefficient T>. As a consequence, the results for the parametric form factor 
can be obtained directly from the spectral form factor without further calculations. One has to 
express the contributions to the form factor, ([16|) and fl36|), in terms of T>\ and T>\ and replace T>* 
by T> 2 . For — 2 one obtains in this way 

K ac (r, x) - K sc (r, 0) « ^ [V{D* 2 — 2iV 1 + 2i V* 2 ] + ^ [{V x V* 2 f 

-4iV l V 2 >(V l - V* 2 ) - AV 1 V 1 - AV* 2 V* 2 + 8V ± V*] 

r 2 t 3 
= — [V{D* 2 - 2i Pi + 2i V* 2 ] + — [V{D* -2iV x + 2i V* 2 f . (40) 

A first point to notice is that the parametric form factor in fl40|) is, in general, not real. This seems 
to be in contrast to the random matrix result (^) which is real. The reason for this lies, however, 
in the correct choice of the unfolding procedure. The definition ( |3T| ) yields only the universal 
form factor in case that the mean density of states d(E) does not depend on the parameter of 
the system. However, in the present case the mean density changes slightly with the parameter 
x of the system, and this leads to a slight shift of the spectrum with x []I9| . As a consequence, 
the argument rj of the two-level correlation function is shifted, and its Fourier transform, the form 
factor, is multiplied by a term of the form e lCT , where c is determined by the shift of the levels. 
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By rewriting ( [40|) up to the considered order in r in the form 

/ t 2 r 3 
K sc {t, x) - K sc (t, 0) w r + — Re [V{D* 2 -2iV x + 2% V* 2 ] + — (Re [V{D* 2 -2iV x + 2% V* 2 ]) 2 

y 4 32 

x exp Im [V{D* 2 -2iV x + 2i D 2 *]j - r , (41) 

one can extract the result that corresponds to a proper unfolding by dropping the exponential, 
and one obtains 

K sc (t, x) - K sc (t, 0) « ^ Re [Z^ - 2z P x + 2* D 2 *] + ^ (Re [2^ - 2z V x + 2i V 2 }) 2 

r 2 t 3 

= |£>i-X? 2 | 2 + |£>i-£> 2 | 4 , (42) 

8 1 1 128 1 1 ' v ; 

where (H) has been used. 

A comparison with (^) shows that it has now the same form as the random matrix result. The 
remaining step is to connect the parameter a with the universal parameter x. Since we consider 
small parameter differences, x is given by x = Aa a v , where a v is the square root of the variance of 
the velocities with respect to the parameter a. In order to evaluate a v we employ a semiclassical 



method |2(J, [21], |22| which expresses it in the form 

*2 = U7rY) =(^vKie)r) • (■13) 



da J I \v-+o 

E 



where dP v is a Lorentzian smoothed density of states, that is weighted by the velocities and expressed 
semiclassically by 



(44) 



7 

Here i] is the width of the Lorentzian and all other quantities are defined as before. 

For small parameter differences we can express the derivative of the diffraction coefficient with 
respect to a by {T> 2 — T>i)/Aa. After inserting ([0]) into ( |4"5D the leading order contribution comes 
from single-diffractive orbits, and we evaluate the double sum in the diagonal approximation with 
the sum rule (0) and amplitudes ( |10|) 

2 i- ^ yMjT-Tj \V 2 -V 1 \ 2 J 2 V T,\ 

*.=j3&y dT E 4ffWl2 | (Aa)2 ^{-m(e)S 

= l ^J dT ^)^A^ eXP |-^)) 
2 |P 2 -P!| 2 

/3 (47rAa) 2 ' 1 J 

For (3 = 2, the universal parameter is given by x = Aaa v = \T> 2 — X>x | / (47r) , and after insertion 
into (f42"D we reproduce the random matrix result ( |39| ) 

K sc (t, x) - K sc (t, 0) ss -2ttVt 2 + 2ttVt 3 . (46) 
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For the GOE ensemble, the parametric density correlation function is given by a triple-integral 
which cannot be expressed in closed form \\L8\ . We consider here only the first-order correction 
for which the GOE-result can be obtained from the asymptotic form of the correlation function 
for long-range correlations in || 23]: K goe (t, x) — K GOE (r, 0) ~ — 2tt 2 x 2 t 2 . Again we find an 
agreement with the semiclassical result K sc (r, x) — K sc (t, 0) ~ — r 2 \D\ — T> 2 \ 2 /4: = —2tc 2 x 2 t 2 . 



6 Discussion 

We have investigated in this article the influence of a point-like scatterer on the spectral statistics 
of quantum systems with chaotic classical limit. It has been shown that the modification of 
the form factor K(r) due to the scatterer can be evaluated systematically in a semiclassical 
approximation. The expansion of the form factor in powers of r corresponds on the semiclassical 
side to an expansion in the number of loops of the diffractive orbits. We have calculated off- 
diagonal contributions to the r 2 - and r 3 -term, but the method can be extended to higher order 
terms. 

The results lead to the conclusion that the delta-perturbation does not modify the form factor. 
Off-diagonal terms from pairs of different diffractive orbits and from pairs of diffractive and peri- 
odic orbits cancel exactly the diagonal terms from diffractive orbits. This requires the existence of 
correlations between different orbits. These correlations arise from pairs of orbits which are very 
close in coordinate space. 

The results provide a support for the random-matrix conjecture. They imply, up to the con- 
sidered order, that the statistics of chaotic systems are invariant under the perturbation by a 
point-like scatterer as is expected from universality. They show also that correlations between 
two energy spectra for different parameter values are universal, provided that the parameter 
difference is small. Furthermore, they indicate indirectly that the spectral statistics of the unper- 
turbed system (and thus also of the perturbed system) are identical with those of random matrix 
theory. The reason for this is that independent results on the invariance of spectral statistics 
under a delta-perturbation are based on the assumption that the unperturbed energy levels and 
wave functions have random matrix distributions E5fl. Since the semiclassical results show this 



invariance for chaotic systems, the combination of both results provides a theoretical indication 
that chaotic systems follow the random matrix hypothesis. 

Finally, the results are a support for the semiclassical method. They show that semiclassical 
approximations are capable to go beyond the leading term in r and are an appropriate tool for 
investigating spectral statistics in the semiclassical limit. 

I would like to thank K. Richter and P. Seba for helpful discussions. After completion of this article 
I learned about work by E. Bogomolny, P. Leboeuf and C. Schmit |?5|| with related semiclassical 
results for the first-order correction. 
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A Stability matrix for the motion in a magnetic field 



The stability matrix M of a trajectory determines infinitesimal orthogonal deviations from the 
final point of the trajectory in terms of the deviations from the initial point of the trajectory 

(£)-*&)■ 

In systems with a magnetic field the momentum has the form p = mv + ^A(r), where q is the 
charge of the particle and A is the vector potential. In this case it is often more convenient to 
consider a matrix M that describes deviations of the velocities instead of those of the momenta. 
The relation between both matrices is given by 

( dTf )=M( dri ) , M = A^MA, A hf =( l °) , (48) 
\mdvf I \mdvi J \ a i,f 1 / 

where a^/ = -(fiij • V)(njj • A(rij)), and fiij is the direction orthogonal to the trajectory at 
the initial and final point of the trajectory, respectively. The matrix M has unit determinant and 
satisfies M12 = M12 and TrM = TrM + M^fai — a/). In cases where the initial and final points 
are identical and the initial and final velocities differ by a small angle e, the traces are identical 
in leading order of e. Since the semiclassical approximations in this article involve only M 12 and 
Tr M — 2, we express all quantities in terms of M instead of M, and we use also the term stability 
matrix for it. 
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